The effect of spin-orbit nuclear charge density corrections due to the anomalous 

magnetic moment on halonuclei 
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In this paper we consider the contribution of the anomalous magnetic moments of protons and 
neutrons to the nuclear charge density. We show that the spin-orbit contribution to the mean-square 
charge radius, which has been neglected in recent nuclear calculations, can be important in light 
halonuclei. We estimate the size of the effect in helium, lithium, and beryllium nuclei. It is found 
that the spin-orbit contribution represents a ~ 2% correction to the charge density at the center 
of the ^Be nucleus. We derive a simple expression for the correction to the mean-square charge 
radius due to the spin-orbit term and find that in light halonuclei it may be larger than the Darwin- 
Foldy term and comparable to finite size corrections. A comparison of experimental and theoretical 
mean-square radii including the spin-orbit contribution is presented. 



I. INTRODUCTION 

It is known that the spin-orbit density of neutrons (and 
protons) contribute to the nuclear charge density as a rel- 
ativistic effect IJ. Consider a reference frame in which 
there is no charge density present. By applying a Lorentz 
transformation, we see that the Q-th component of the 4- 
current vector (the charge density) develops a non-zero 
value if the vector current is non-zero. The vector cur- 
rent exists in neutrons due to the anomalous magnetic 
moment. 

The spin-orbit density effect of neutrons has been 
shown to contribute to the change in nuclear charge ra- 
dius (and hence scattering cross-sections) between ^°Ca 
and ^^^Ca [l| (see also for more recent relativistic cal- 
culations). On the other hand the spin-orbit contribu- 
tion was seen to have a negligible effect on the differen- 
tial cross-section of ^"^Pb J^, and by extension other 
heavy nuclei. 

The relativistic correction known as the Darwin-Foldy 
term is by convention a recoil correction. In a recent 
review paper Q , it was recommended that this correction 
be treated as part of the charge radius, due to the fact 
that it appears in the charge density of the proton. This 
will be the approach we will adopt in our calculations. 

Interest in the neutron contribution to nuclear charge 
densities has been renewed by the discovery of halonu- 
clei. Comparison between theory and high-precision ex- 
periment has been performed in isotopes of He 0, Q and 
Li [ol-fisj. Relativistic corrections and finite size effects 
are consistently taken into account in these theoretical 
calculations, while the spin-orbit effect is never explicitly 
considered. 

In this paper we calculate the spin-orbit contribution 
to nucleon charge density. We derive a simple, general 
expression for the contribution of the spin-orbit density 
to the mean-square charge radius. In light halonuclei this 
contribution may be comparable to the finite size effects. 



and may be larger than other relativistic corrections that 
are included such as the Darwin-Foldy term. In addition, 
we use the oscillator model to obtain numerical estimates 
of the spin-orbit charge density in ^Be. 

We present the calculated values of the mean-square 
charge radii, including the spin-orbit contribution, in Ta- 
ble IVl 



II. SPIN-ORBIT CHARGE DENSITY IN THE 
NONRELATIVISTIC LIMIT 

We begin with the expression for the 4-current derived 
from quantum field theory considerations [Tg'I (relativis- 
tic units h — c — 1 are used throughout, unless otherwise 
stated). 



q^F„ 



where M is the mass of the nucleon, and F„ 



7^ 2M \ 

(1) 

are elec- 
tromagnetic form factors, = 4M^ — q'^ is the momen- 
tum operator, q is the momentum transfer operator, and 
^fii^ _ _ j^jf^'^ is the antisymmetric tensor. For 

nucleons at the Fermi surface jjj ~ ^ 0.08, there- 
fore in this paper we perform all calculations in the non- 

2 

relativistic approximation, accurate to order -p^. Here 
Pp refers to the nuclear Fermi momentum. In the non- 
relativistic limit the lower component of the Dirac wave- 
$ 

X 



function ip — 



is given simply by 
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The charge density can be calculated from the 0-ift- 

is 



component of the 4-current ([T]), which to order 
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The form factors and Fm are, in general, functions 
of the momentum transfer squared . However, for our 
purposes it suffices to consider only Fe(0), the electric 
charge of the nucleon, and F,„(0) — Fe(0) = /i, the anoma- 
lous magnetic moment — —1.91 and /ip — 2.79 — 1 = 
1.79 for neutrons and protons, respectively). 

Substituting the nonrelativistic wavefunction into ^ 
gives the spin-orbit contribution to the charge density, to 



order 



4M2 



(4) 



Here there is an implied summation over a,/3 = 1,2 — 
the upper and lower components of the spinor Note 
that the Darwin- Foldy term [F^ /8M'^) can also be de- 
rived from (|3]). 

Calculation of the spin-orbit charge density is de- 
scribed in the Appendix. 



where 



-(/ + 1), J=l + ^ 



(7) 



We defer proof to the Appendix. Interestingly, {r'^)i]} is 
independent of the radial wavefunction /(r). 

To compare calculated point-proton mean-square 
charge radii with experiment, we introduce nuclear size 
and relativistic corrections Q using the equation 

N 

{r^)cH - (r2)pp + {RD + -{Rl) + ^ + {r^)so (8) 



where 



III. CHARGE DENSITY IN SPHERICALLY 
SYMMETRIC POTENTIALS 

Assuming a spherically symmetric nuclear potential, 
we can write the wavefunction in the shell-model as 
^ = f{''^)^jim, where ft refers to the spherical har- 
monic spinors (jA.2[) . From Q we find that the anoma- 
lous charge density produced by s-waves is identically 
zero. This can be readily seen by noting that for each 
■Si/2 projection, the wavefunction only has either an up- 
per or lower component, which means that summation 
over a and (3 selects only the (Tz.aa component (in the 
standard representation). Since the wavefunction com- 
ponents in this case are real, V<i>* = V$q, and therefore 
Ta/J • V$l, X = 0. 

Consider what happens to the total local charge den- 
sity when we sum over a spin-orbit doublet of nucleons. 
In the LS-coupling scheme, the (Tai3 term in Eq. re- 
sults in the cancellation of p for pairs of antiparallel spins. 
This implies that the sum over any spin-orbit doublet is 
zero. On the other hand in the JJ-coupling scheme there 
is no straightforward cancellation, and as a result the sum 
of the charge density contribution over a closed subshell 
is not zero. However, the sum over a spin-orbit doublet, 
for example IP3/2 and lpi/2, is zero if the radial wave- 
functions are equal (any difference is of higher order in 

M'- 

For comparison with experiment, one wishes to know 
the contribution of the spin-orbit charge density to the 
mean-square charge radius. We define this for a single 
nucleon as 



2\(1) _ 



Using this definition and Q, we find 

,2\(1) _ M 



(k + I) 



(5) 



(6) 



and i runs over all nucleons (but there is zero con- 
tribution from a sum over spin-orbit doublets). Here 
{Rp = 0.769(12) fm^ and (Rl) = -0.1161(22) fm^ [l^ 
are the mean-square charge radii of the proton and the 
neutron respectively, and the term = 0.033 fm^ is 
known as the Darwin- Foldy term [(|]. The effect we are 



considering simply adds the quantity 
tion used by previous works, e.g. [l8| . 
from normalization. 



(r )so to the equa- 
The factor ^ arises 



IV. SPIN-ORBIT CHARGE DENSITY IN 
BERYLLIUM- 7 

As was shown in the previous section, {r'^)so does not 
depend on the radial wavefunctions. However, one may 
need to know the value of the charge density p{r), which 
depends on the radial wavefunctions. To estimate the size 
of the effect, we will use the harmonic oscillator potential 



y(r) = 



(10) 



as an approximation to the nuclear potential. Using the 
known solutions of the quantum harmonic oscillator [l9| , 
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the form of the spin-orbit charge density pjim (jl]) for the 
Z = 1, j = 3/2 states can be written as follows 



= = -^A\r){l - 2a{x' + i/)) (13) 
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2^2 2 ^ 2 
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A^(r)(3-2a(x"+y^ + 4z^)) 

(14) 



3 



TABLE I. Comparison of '^Be experimental and theoretical 
charge radii in fm^. 

{r^)exp [201 (r')NCSM [15J (r^)GFMC [8] {r^)so 
^Be 6.294(0) 5.83 6.67 0.02 



where A'^{r) = ^l^p-^""'^ and a = ^ parametrizes 

the potential. Note that the sum of these charge densities 
is spherically symmetric: 

T.Piim = (15) 

m 

a can be obtained through the fitting of the model to 
known experimental measurements. 

To estimate the size of the spin-orbit charge density 
in ^Be, we use the shell model, where the neutrons and 
protons fill nuclear states independently. ''Be has a single 
unpaired neutron in the lp3/2 state. The charge density 
generated by this neutron is given by (|13p. By compari- 
son, in our model the usual charge density contribution 
from the four protons, taking into account only highest 
order terms, is (e = 1 in our units) 

pp,totai = y^^i^^ 

i 

= ^(x^+.^ + ^)e--^ (16) 

The relative contribution of the spin-orbit charge density 
to the total nuclear charge density is the greatest at the 
center of the nucleus, where it represents a 2% cor- 
rection. For this calculation, we used uj = ^'^j^ij^^ as our 
fitting parameter for the harmonic oscillator potential; 
this corresponds to a = 0.25 fm~^. 

The effect of the spin-orbit charge density on the mean- 
square charge radius is given by (jH]) . For the neutron and 
two protons in the lp3/2 state it is 

{r^)so = ^-^(Ai„ + 2/ip) 

= -0.021 + 0.039 = 0.02 fm^ 

This is a maximal estimate; configuration mixing reduces 
the size of the effect. Table U presents the experimental 
and theoretical charge radii of ^Be. It shows that the size 
of {r'^)so can be important, although it is smaller than the 
spread between different theoretical data points. 

V. HELIUM HALO NUCLEI 

There has been a considerable increase of interest in 
the properties of halonuclei, as evidence d by the number 
of recent papers written on the topic [l8l l20l - |22| . In *He, 
the neutrons are assumed to completely fill all 4 states 



TABLE II. Comparison of ' He experimental and theoretical 
charge radii in fm^. 





{r2)GFMC M 




''He 4.277(0) [23| 4.04 


4.15 


-0.08 


*He 3.722(0) [21J 3.89 


3.65 


-0.17 



of the lp3/2 orbital. The calculated change in the mean- 
squared charge radius due to this effect is then 

(^')so = E / ^^Plirnd'r = -0.17 fm^ . (17) 

m 

This shows that the contribution of a closed subshell can 
be significant. It is interesting that while the charge den- 
sity for different projections appear to be different, the 
correction to the charge radius is independent of projec- 
tion for any fixed k. 

The spin-orbit term contributes a decrease in the 
charge radius from ^He to *He, in addition to that pre- 
dicted by the point-proton model. The magnitude of this 
effect can be much larger than that of the Darwin- Foldy 
term in the relation ([S]), which justifies the inclusion of 
this spin-orbit effect. Table |lT] shows current available 
theoretical (no-core sheh model (NCSM) and Green's 
function Monte Carlo (GFMC) [sl) and experimental 
data, and the estimated size of the spin-orbit effect. We 
note that the calculated spin-orbit contribution is a max- 
imal estimate since configuration mixing may reduce the 
effect. 



VI. LITHIUM NUCLEI 

It is easy to extend our results from ^He to ^Li be- 
cause in the shell model the neutrons fill the same states. 
The unpaired lp3/2 proton in the lithium nucleus then 
provides the following correction to the charge radius 

{r^)so,p = 0.03 fm^ . (18) 

Experimental results [31 for ^Li and ^^Li show that ^^Li 
has the bigger charge radius. For neutrons, our calcu- 
lated spin-orbit contribution is negative for the ^Li nu- 
cleus. Because the neutrons form a closed p-shell in "'^^Li, 
their spin-orbit contribution in this case is zero. Table 
mil also shows that the size of the spread in theoretical 
data points is much larger than the spin-orbit correction, 
therefore it is not possible at this moment to comment on 
the effect that the inclusion of this term has on the accu- 
racy of theoretical results. However, the size of this effect 
can be larger than existing corrections that are taken into 
consideration, and therefore should be accounted for as 
well. 
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TABLE III. Comparison of ^'^'^'i^Li experimental and theoretical charge radii in fm^. 





) exp 


(r2)GFMC \g,M 


(r')svMC [11,22] 


(r2)DCM [13] 


(r")NCSM [44,_15] 




«Li 


6.336(7) [24| 


6.40 




7.19 


5.08 


-0.00 


^Li 


5.726(4) [24] 


5.71 


5.80 


6.46 


4.83 


-0.03 


»Li 


5.286(8) [25] 


4.98 


5.36 




4.38 


-0.06 


»Li 


4.916(6) \18} 


5.37 


4.98 


6.43 


4.36 


-0.09 


"Li 


6.087(9) fl8] 




6.40/5.ir 


7.36 


4.44 


0.03 



" The second value assumes a frozen Li core. 



TABLE IV. Comparison of experimental charge radii with 
current GFMC calculations (the value for ^^Li is from 
NCSM). The second column represents current calculations 
with the spin-orbit effect we calculated taken into account. 
All values are in fm^. 





('■^)gfmc 




'} exp 


''He 


4.15 [8] 


4.07 


4.277(0) [23| 


«He 


3.65 [8] 


3.48 


3.722(0) [21] 


«Li 


6.40 \9,W} 


6.40 


6.336(7) \2£ 


^Li 


5.71 [9,20] 


5.68 


5.726(4) \2£ 


«Li 


4.98 [9,20] 


4.92 


5.286(8) [25| 


»Li 


5.37 [9, 10} 


5.28 


4.916(6) [18| 


"Li 


4.44 [14,15] 


4.47 


6.087(9) [18] 


^Be 


6.67 [8] 


6.69 


6.294(0) [20| 



VII. CONCLUSION 

We have calculated the spin-orbit contribution to the 
charge density of nucleons. In ^Be, the spin-orbit charge 
density contributes approximately 2% to the total charge 
density at the centre of the nucleus. The effect of the 
spin-orbit density on the mean-square charge radius has 
been shown to be independent of the radial wavefunc- 
tions, and hence the form of the (spherical) potential. A 
summary of the results of our calculations for the mean- 
square charge radius can be found in Table IIVI 

In halonuclei such as ®He and ^Li, our results show 
that this effect presents a considerable contribution to the 
charge radius, sometimes larger than the Darwin-Foldy 
term and comparable to the finite size effects that are 
already taken into account. It is worth noting, however, 
that our calculations do not include configuration mixing 
which may reduce the size of the effect. We suggest that 
the spin-orbit contribution should be calculated in the 
same scheme as the point-proton radii. 



Appendix 

Here we prove that the contribution of the spin-orbit 
charge density to the mean-square charge radius of a nu- 
cleus due to a nucleon in a spherically symmetric nuclear 



potential is given by Eq. ^ . We start with the nonrela- 
tivistic expression (jH): 



P 



2i(T^l3- [(VO X (V$,3)] 



where there is an implied summation over a, /3. This 
equation can be derived from ([3]) by substituting the non- 
relativistic form of the lower component of the wavefunc- 
tion (l2|). Assuming a shell model allows us to write 



* = /(r) %i™(e,0) 



jlm 



/ m-t-i, i 



1 1 iXu 



(A.l) 
(A.2) 



where /(r) is the (real) radial wavefunction, fljimiO, (j)) 
is the spherical harmonic spinor, and the C are Clebsch- 
Gordan coefficients relating the orbital angular momen- 
tum I and the spin to the total angular momentum j. 
The gradient of this wavefunction is then 

V$ = r2 /,. f + / Vf7 . 

where = df/dr and we have dropped the subscripts 
from ft. Substituting into (H]) gives 



4Af2 



2i (T. 



—nUr X vr!«) - ^(r X vr!„)tr!« 

r r 



f ynl X 



(A.3) 



We split this into two parts which we treat separately: 
p = pi + p2, the first line and second lines of (IA.3|) . 
respectively. Using the relationship 



l(f> = (r X p)$ = ~ifr X Vfl 
the first line of (jA.31) can be written as 

r.^ fr 



Pi 



2/1 
4Af2 


r 


2fi 




4M2 


r 





I' 



(A.4) 
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where = j(i + 1) and similar for / and s. The contri- 
bution of p to the mean-square charge radius is given by 
([5]). Therefore the contribution of pi is 



K'P ' J 2 dr 



(A.5) 



where integration by parts was used in the last step and 
f{r)ft{9, 0) is assumed to be normalised to unity. 

The second line of (|A.3p poses more of a challenge. 
Using the identity 

r 0& r sm 9 ocp 

where (f , 6, <f)) form an right-handed orthogonal basis, 
we write the contribution of p2 to the mean-square charge 
radius as 



2Af2 
2M2 



de 



de 



The Pauli matrices can be projected onto the spherical 
basis vectors: cr = + agO + a^cf) with 



cos t/ sm f e 

sin 9 e"^ — cos ( 



-20 



= 



— sm tl cos y e 



COS tie 

- 



sin 9 







(A.6a) 
(A.6b) 
(A.6c) 



Substituting Gr and summing over the indices a and /?, 
we can write 



2A'P 



COS 9 



sm V e 



sm a e 



— cost 



I d9 d<j) del) 89 

yd9 d(f> ~ d(j) 89 j 
'dVLldVL2 dill 8^2^ 



d9( 



2A'P 



{h +h+h + h 



(A.7) 



where f2i and VL2 refer to the upper and lower compo- 
nents of VL, respectively, and the four integrals will be 



considered separately. The spherical harmonic functions 
in the definition of $7 (|A.2p can be expressed in terms of 
associated Legendre polynomials, P™: 



which allows us to write Vl as 

^^/c,Pr?(cos.)e^(™-^)^\ 
V C2Pr^ (cos 0) e^(™+^)^ ; 

where Ci and C2 are real coefficients. Using this definition 
of f2 and the transformation u = cos 0, we can express 
the integrals (|A.7p as 

h = 2.27r c?(m - -) / mP""^ (u)'P""^ 

h = 27rciC2 (^(m + i)P""^ (u)'P;"+' (^.) 

h = h 

h = ~2.27r c2(m + ^) / mP;"+^ (u)'P;"+^ 

where Pl^{u)' refer to derivatives of the associated Leg- 
endre polynomials. 

Ii and /4 can be solved using integration by parts and 
the orthogonality relation 



2(lA-mV 

Pr{u)Pnu)du = ^^_i_4, . (A.9) 
1 [2l + l){l-m)\ 



I2 and /a are equal to each other and can be solved by 
transforming the derivatives of the polynomials using the 
identities 

(W^ - l)P™-5(u)' = Vl -U2p™+5(y) 

+ (m- V2)MPr~^(u) 
(u^ - 1)P"+^(m)' = -{l + vi+^/2){l - m + 1/2) 



• VT^^P" ^ (u) - (m + V2)wPr ' (u) 



and then using the orthogonality relations (jA.9|) . The 
results are 



= 1/2) {ci2^,,,f 



h - i iCf;:_^i i _,va-"^+V2)(^+m+v2) 

2 ' 2 2 ' 2 ' 2 2 

i4 = (m + V2)(q'" )^. 
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Substituting into (IA.7I) and evaluating the Clebsch- 
Gordan coefficients we obtain 



(A.IO) 



Adding this contribution to (jA.5[) gives the total change 
in mean-square charge radius due to a single nucleon 

The contribution of a filled subshell to the normalised 
mean-square charge radius is therefore 



1 



7 E <^)^ 



1) _ 2i + 1 



subshell 



Z A'P 



[n + l) 



in agreement with [ij. For a spin-orbit doublet, on the 
other hand, we find 



SO doublet 
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